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Abstract 

We address the flavour composition along the border between the hadronic 
and the quark-gluon plasma phases of QCD. The ratio of strange to up 
and down antiquarks (A s ) produced in particle and nuclear collisions, is 
found to increase in collisions with initially reached energy density (ej) 
up to e cr it ~ 1 GeV//m 3 . Above this value it decreases approximately 
linearly and reaches its asymptotic value at zero baryon chemical poten- 
tial (/is). We demonstrate that X s in nuclear collisions is approaching 
its asymptotic value at e« ~ 8-9 GeV//m 3 , corresponding to ~ 3- 
8 TeV per nucleon+nucleon pair which will be reached at the LHC. 
After correcting for the difference in the chemical potentials of vari- 
ous colliding systems, A s universally saturates across the QCD phase 
boundary, following the temperature. Recent experimental puzzles as 
the increase in the K/ir ratio in Pb+Pb collisions at 40 GeV per nu- 
cleon, its different behaviour at midrapidity, the decrease of the double 
ratio of i*T/7r(A+A/p+p) in nucleus nucleus over p+p collisions with 
increasing y^s, and the increase of \ s in p+A over p+p collisions at the 
same i/s, are naturally explained. We study the approach of thermo- 
dynamic observables at hb = to the transition point and extract an 
estimate of the critical temperature. 



We thank the Schweizerische Nationalfonds for his support. 



1 Introduction 



One outstanding prediction of the theory of strong interaction (e.g. 
which applies to the evolution of the early universe, is the phase transition from 
confined hadrons to a deconfined phase of their constituents, the quarks and 
gluons (the so called quark gluon plasma =QGP). An experimental program 
which started approximately in the eighties and continues in many accelerators 
as: CERN SPS, BNL RHIC and others, has been dedicated to the experimental 
verification of this transition |§||0@fl]' ^he ma i n process is nuclear reactions 
at high energy to achieve a thermalized state with temperature exceeding the 
critical temperature for the QCD phase transition. One of the predicted sig- 
natures of this change of phase, is an enhancement of strange particles [fLOf 



The main idea was a) the lower threshold for the production of s~s quarks in 
the QGP through e.g. gg — > ss as compared to the higher threshold for pro- 
duction of strange hadrons in a hadronic reaction for example pp — > pAK + 
and b) the similarity of the mass of the strange quark with the critical tem- 
perature T ~ 200MeV, which allows for the equilibration of strange quarks in 
the QGP. 

Argument a) holds for all quark flavours, leading to the conclusion that all 
hadron multiplicities are expected to be increased when produced out of a 
hadronizing QGP, as compared to hadron production out of a collision which 
does not pass through the QCD phase transition. The enhancement can be 



related in a simple way to phase space [ 'I I. ] and is flavour dependent . 



However the argument b) is true only for the quarks which have mass less or 
of the order of the critical temperature of approximately 200 MeV PJ . 
Therefore next to the up and down quarks, strange quarks are expected to play 
a crucial role in identifying the QCD phase transition. In particular, hadrons 
with u, u, d, d, s and s quarks are expected to reflect the critical behaviour of 
a locally equilibrated phase made up of u, u, d, d, s,s quarks and gluons. 
Heavier flavours like charm and beauty can be affected by the QCD transition 
in several ways discussed in the literature P~2"|,|8|,[T3| . They give important in- 



sights e.g. by the experimental measurement of the dissociation temperature 



of cc or bb states which can be at or above the critical one [12], 
Therefore, since dissociation of quarkonia may be due to overcritical energy 
densities, it appears that the critical parameters of the transition can be ex- 
tracted only through other observables which are sensitive to the transition at 
the critical energy density and approximately frozen with fast hadronization. 
For example fluctuations of several parameters and multiplicities of hadrons 
which include up, down and strange quarks. 

In this paper we concentrate on the behaviour of the strange flavour across 
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the QCD phase boundary. The main point is to interpret the global features 
of the data (A+A, pp, e + e~ collisions at -^5=2-1800 GeV) within a strictly 
thermodynamic approximation, and not to reproduce the data in detail using 
models outside the validity of thermodynamics. 

Clearly, no perfect global equilibration is reached in these systems which more- 
over are a mix of several thermodynamic reservoirs e.g. along the rapidity axis. 
It appears however that the hadrons in the final state of the colliding systems 
studied, do not rule out a thermodynamic description. This is especially the 
case when a centrality selection has been imposed. 



Preliminary results of part of this work have been shown in |L4| 



In section 2 we extract thermodynamic parameters (T, /xb,/z s ) for fixed target 
Au+Au collisions at 2 and 4 GeV per nucleon and for fixed target Pb+Pb col- 
lisions at 40 GeV. In section 3 we estimate the initial energy density reached 
in the collisions. In section 4 we examine the strangeness suppression factor 
at non zero chemical potentials as a function of the initial energy density, and 
find its asymptotic value using nuclear collisions. In section 5 we extract the 
temperature at zero chemical potentials and draw the 'strange border of the 
QCD phases'. We discuss how the critical parameters for the transition are 
extracted and elaborate on the explanation of several recent experimental puz- 
zles. 



2 Thermodynamic description of nuclear col- 
lisions at 2, 4 and 40 A GeV 

We compare the ratios of experimentally measured hadron yields in nuclear 
collisions with the prediction of a grand canonical ensemble of non interacting 
free hadron resonances. We consider the pseudoscalar and vector u,d,s meson 
nonets as well as the spin 1/2 baryon octet and spin 3/2 decuplet and their an- 



tiparticles as well as the /o(400 — 1200) or cr, interpreted as scalar glueball |LS 
We impose exact conservation of strangeness and correct for the change of final 
observed hadron yields due to decays of higher lying resonances (the so called 
'feeding'). Further details of this model can be found in |4j]. We extract then 
the thermodynamic parameters: temperature, baryochemical potential (/4b) 
and strangeness chemical potential (fis) from the model prediction which de- 
scibes the data best. We discuss in the following the production of strangeness 



(25) 
{u+d) 

a measure of the strangeness suppression factor defined and used in the liter- 



relativ to non strange particles extracting the quantity: X s = wn i cn ^ s 
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ature as A s 



2(s+s) 



(u+u+d+d) 

Similar analysis extracting thermal parameters from data, can be found in the 
literature e.g. 



1710. 8119,16, 20, 21, 22 



24 



33| . For a recent review see 
A main new idea introduced in [11.3J] and used also in this paper as a tool, 
is the extrapolation of all thermodynamic states to equivalent states at zero 
chemical potentials e.g. along isentropic paths. 



Throughout the paper, all parameters discussed concern the state of hadrons 
at their chemical freeze-out, that is the time after which the hadron yields 
do not change anymore through inelastic collisions. The only exception is the 
initial energy density which refers to the early state of the system. 



2.1 Thermodynamic description of Au+Au collisions at 
2 A GeV 

We use 3 ratios measured in fixed target Au+Au collisions at 2 GeV per nucleon 
beam energy shown in table p] from references |25|j26]| and impose strangeness 
conservation. The predicted and the experimental ratios are shown in table |l| 
The K/tt and the n/p ratios are measured at midrapidity, while the A/ K® ratio 
nearby p5| . To account for the different phase space we add a 10% systematic 
error to the ratio. 

The resulting x JDOF is (1.41/1) (CL 23%). After defining the (T, fi b , fj, s ) 
values describing the particle ratios produced in central Au+Au collisions at 
^=2.3 GeV at the chemical freeze-out we extrapolate to the T at zero fugac- 
ities (table §) along an isentropic path. 



ratio 


model 


data 


K/tt 


0.0110 


0.00929 ± 0.00257 


Tv/p 


0.268 


0.266 ± 0.0770 


A/K° s 


2.770 


3.222 ± 2.274 



Table 1: Au+Au at y/s= 2.3 GeV. 

Predicted versus experimental particle ratios for the best fit of our model. 



2.2 Thermodynamic description of Au+Au collisions at 
4 A GeV 

We use 4 ratios measured in fixed target Au+Au collisions at 4 GeV per nucleon 
beam energy shown in table [5] from references |25|.[26| and impose strangeness 
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Ps 


T 


Xs 


Ps 


T(eq,p s ) 


X s (eq,p s ) 


GeV 


GeV 


GeV 




1/fm 3 


GeV 




0.714 


0.0918 


0.0510 


0.0321 


0.0435 


0.072 


0.0388 






± 0.006 


± 0.0327 




+0.032 -0.013 


+0.116 -0.0249 



Table 2: Au+Au at ^= 2.3 GeV. 

Thermodynamic parameters for the best fit and temperatures and A s extra- 
polated to zero fugacities. 



conservation. The predicted and the experimental ratios are shown in table [5| 
The K/tt, tt/p and K + /K~ ratios are measured at midrapidity, while the 
A/if° ratio nearby [p5| . To account for the different phase space we add a 10% 
systematic error to the A/K® ratio. 

The resulting x 2 /DOF is (1.62/2) (CL ~ 40%). After defining the (T, (i b , 
/i s ) values describing the particle ratios produced in central Au+Au collisions 
at t/s=3.3 GeV at the chemical freeze-out we extrapolate to the T at zero 
fugacities (table |) along an isentropic path. 



ratio 



model data 



K/w 0.0396 0.0422 ± 0.0111 

K+/K~ 12.437 12.316 + 0.699 

K/K° s 2.342 2.750 ± 0.342 

n/p 0.473 0.435 ± 0.114 



Table 3: Au+Au at ^s= 3.3 GeV. 

Predicted versus experimental particle ratios for the best fit of our model. 



Pb 


/', 


T 


A, 


Ps 


T(eq,p s ) 


\s(eq,p s ) 


GeV 


GeV 


GeV 




1/fm 3 


GeV 




0.657 


0.0923 


0.073 


0.139 


0.163 


0.097 


0.125 






± 0.005 


+0.047 -0.040 




+0.016 -0.009 


+0.0710 -0.0359 



Table 4: Au+Au at ^/s= 3.3 GeV. 

Thermodynamic parameters for the best fit and temperatures and A s extra- 
polated to zero fugacities. 
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2.3 Thermodynamic description of Pb+Pb collisions at 
40 A GeV 

We use 5 ratios measured in Pb+Pb collisions at 40 GeV per nucleon beam 
energy shown in table |5| from references and impose strangeness con- 

servation. The predicted and the experimental ratios are shown in table 
The ratios A/A, 5/H are measured at midrapidity, while the other ratios are 
given in full phase space acceptance. The (B — B) is taken equal to the total 
number of participant nucleons. 

To account for the different phase space and the use of the assumption (B—B) = 
^participant we add a 10% systematic error to the 7r/(B — B), A/A and S/H 
ratios. 

The resulting x 2 /DOF is (6.68/3) (CL ~ 10%). After defining the (T, fx b , 
Us) values describing the particle ratios produced in central Au+Au collisions 
at ^=8.76 GeV at the chemical freeze-out we extrapolate to the T at zero 
fugacities (table along an isentropic path. 



ratio 


model 


data 


K/ir 


0.128 


0.125 ± 0.0072 


K+/K- 


3.0355 


3.163 ± 0.232 


n/(B-B) 


1.122 


0.852 ± 0.125 


A/A 


0.0216 


0.0230 ± 0.00251 




5.146E-2 


8.00E-2 ± 2.625E-2 



Table 5: Pb+Pb at y/s= 8.76 GeV. 

Predicted versus experimental particle ratios for the best fit of our model. 



Mb 




T 


A s 


Ps 


T(eq,p s ) 


\ s (eq,p s ) 


GeV 


GeV 


GeV 




1/fm 3 


GeV 




0.405 


0.090 


0.150 


0.660 


2.466 


0.164 


0.407 






+0.019 -0.039 


+0.138 


-0.326 


+0.024 -0.049 


+0.068 -0.199 



Table 6: Pb+Pb at y/s= 8.76 GeV. 

Thermodynamic parameters for the best fit and temperatures and A s extra- 
polated to zero fugacities. 



5 



The correlation between the extracted thermodynamic parameters tempera- 
ture (T), strangeness suppression factor (A s ) and the chemical potentials of all 
systems together with results from reference H, is shown in the figures [TJ, ^ 
i § § 



3 Initial energy density estimation 

We estimate the initial energy density for the collision systems studied in 
section 2 by taking the nuclear energy density of two overlapping nuclei 2 
times the 7 factor of the colliding particles in the center of mass minus one: 



e 7 = 2 e A (7 - 1) (1) 

with 7 = (y/s/2)/m nuc i eon , and 6^=0.138 GeV/fm 3 is the normal nuclear 
matter density. The value in equation [I| multiplied by the stopping power gives 
an estimate of the initial energy density available for heating. This formula is 
better suited for low energy data for example AGS where the applicability of 
the Bjorken formula for the energy density calculation is questionable. For 
a discussion of the systematic error of ~ 30-50% on the initial energy density 
calculation for nucleus nucleus and particle collisions respectively, see reference 
P,|29[|. The resulting initial energy densities are: 



e in (Au+Au v/i=2.3 GeV) = 0.12 GeV//m 3 , 
e in (Au+Au v^=3.3 GeV) = 0.21 GeV//m 3 , 
e in (Pb+Pb v/i=8.8 GeV) = 1.01 GeV//m 3 , 



4 The energy density dependence of X s at fi- 
nite (i 

Two general comments are important for the understanding of the behaviour 
of strangeness in particle and nuclear collisions and our later discussion: Since 
collisions of protons and collisions of nuclei at the same energy per nucleon 
do reach different initial energy densities, it is proper to compare hadronic 
observables as a function of the initial energy density instead of yfs. 
It is also a fact that in all investigated colliding systems with a initial nonzero 
net baryon number, the baryon and strangeness chemical potentials are nonzero 
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and different. Next to the initial baryon number, also the energy of the colli- 
sions and hence the stopping influences the baryochemical potential of the final 
hadrons produced. Therefore the final state of the same nuclei, when colliding 
at different energy, will be described by a different baryochemical potential. 
This leads to the conclusions that a) the behaviour of strangeness as a QGP 
signature can be discussed while comparing systems with the same baryon 
and strangeness chemical potential and b) the initial energy density is a better 
'critical' parameter than y^s, against which different colliding systems can be 
compared. We take these as starting points and study their consequences in 
the following discussion. 



When comparing the strange to non strange particle ratio in nucleus nucleus 
collisions with p+p collisions at the same energy an enhancement is seen, which 
increases with decreasing energy ||3(J. This observation is quantified in figure |6] 
taken from reference [p8|1 , which shows the ratio of ( A + A + K ) /ir as a function 
of the variable F which is a function of the \fs of the collision. The line shows a 
model prediction from reference |31]| assuming that the phase transition occurs 
at the vicinity of the maximum. (For recent literature on the interpretation of 
data on strangeness see e.g. [p^] .) 

In figure |7| the strangeness suppression factor X s is shown as a function of the 
initial energy density. The open stars show nucleus nucleus collision results 
from the present analysis and the ones from reference Q , while the closed stars 
show results from pp and e + e~ collisions from reference M. The lines shown 
are linear fits to the data points to illustrate the tendency of the data. 
The X s points of all investigated nuclear collisions, define the X s behaviour at 
non zero and varying chemical potentials (the two lines (a) and (fl) of the 
upper half of the triangle shown). The horizontal line (7) going through the 
pp and e + e~ data defines the A s value at the zero chemical potentials in the £j 
region shown. 

The \ s factor grows reaching a maximum around ~ 1 GeV/fm 3 and then 
decreases continuously and almost linearly towards an asymptotic value which 
is defined by the curve with zero chemical potentials of the pp and e + e~ data. 
Therefore it follows that no saturation of the X s value is reached at energy 
density between 1 GeV/fm 3 and 8 GeV/fm 3 in nucleus nucleus collisions, as 
the one indicated in figure 10. 



We will argue that the phenomenon of enhancement of the X s value in nucleus 
collisions as compared to the X s values in elementary particle collisions, quan- 
tified in figure |^ by the area of the triangle, is mainly due to the difference in 
the baryochemical potential. 

With increasing energy, the central rapidity region in nuclear collisions becomes 
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increasingly net baryon free approaching a state of a nucleus+antinucleus col- 
lision. 

This tendency is shown in figure |7| by the decrease of the line (a) with in- 
creasing €{. The statistical significance of the present data does not allow a 
thorough investigation of the exact shape of the decrease. In particular in 
order to find deviations from a linear behaviour, more data are needed. 
The nonzero potential A s line (a) crosses the zero potential X s horizontal line 
(7) at an initial energy density of ~ 8-9 GeV/ fm 3 (depending on the way we 
extrapolate e.g. using a linear, exponential or a polynomial distribution). The 
linear fit crosses at 8.25 GeV/ fm 3 . 

This indicates that the limiting value of X s (and an almost net baryonfree 
midrapidity region) will be achieved at nuclear collisions reaching an initial 
energy density of ~ 8-9 GeV//m 3 .0 Assuming that scales as the logarithm 
of y/s, we find that the a/s needed to achieve the limiting value of A s using 
nucleus nucleus collisions is approximately 3-8 TeV. This energy density could 
therefore be achieved by the Large Hadron Collider at CERN. 

5 The energy density dependence of X s at zero 

In the following we discuss the dependence of the X s factor after extrap- 
olating all thermodynamic states with non zero baryochemical potential to 
zero, along an isentropic path. The result is shown in figure |] in logarithmic 
representation (to show the point of smallest A s ), and in figure [|in linear rep- 
resentation. All X s values show a universal behaviour increasing from below 
until 6j ~ 1 GeV/fm 3 and then saturating. 

Several comments can be made: 

1) The figure clearly demonstrates that the peak of X s and of the K/tt ratio 
at 40 A GeV Pb+Pb collisions and their drop towards higher energies, seen in 
figures ^| and [F] is due to the nonzero baryochemical potential in the collisions. 
The X s factor is the same for Pb+Pb collisions at 40 and at 158 A GeV in 
figure |^ (4th and 8th points from the left). 

Therefore, the so called 'strangeness suppression' phenomenon p8[|33], that is, 
the decrease of the strange to non strange particle ratios (e.g. K/tt) from 40 
GeV Pb+Pb towards 158 GeV Pb+Pb collisions is explained by the different 

2 The net baryonfree midrapidity region in a particle+particle collision and its equivalence 
with particle+antiparticle collision, can probably be achieved only in the limit of infinite 
energy. We address here the question at which energy this asymptotic value is significantly 
approached within the errors. 
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baryochemical potentials of these systems. 

This interpretation is supported by the difference seen in the i/i dependence 
of the K + /tt + and the K~ /ir~ ratio as a function of [55]. In particular it is 
seen that the peak at 40 A GeV Pb+Pb colisions, appears only in the K + /tt + 
ratio (K + is 'forced' by the high fis to be abudantly produced in association 
with A). 



2) The increase of the double ratio K/tt (A+A/p+p) with decreasing -^i dis- 
cussed e.g. in the review talk |54j] is a natural consequence arising from the 
above ideas.| 

From this comparison the apparent strangeness enhancement is increasing to- 
wards lower energies. However the comparison leading to this conclusion is 
not taking into account the varying characteristic parameters. 



3) The K/tt ratio is investigated in reference |33[ and shows the above men- 
tioned strangeness suppression with y/s, when constructed using yields of par- 
ticles in the full phase space acceptance, while it remains constant at midra- 
pidity. 

This can be explained by the fact that the baryochemical potential reaches its 
lowest value at the midrapidity of nucleus nucleus collisions and the highest one 



in the forward and backward rapidity regions (multiple reservoirs) e.g. [p5|,p€ 
Therefore, when considering K/tt at midrapidity the bias from nonzero is 
minimized as compared to full acceptance yields. 



4) The enhancement seen in strange to non strange particle yields in central 
p+A collisions when compared to p+p at the same energy [|37|], is explained 



in the same way, since central p+A collisions reach a higher /i and a higher 
initial energy density as p+p collisions at the same a/s. 



5) The strangeness enhancement as defined in the literature (double ratio of 
X s in A+A over p+p collisions at the same -y/i) and as illustrated in figure || 
finds the same explanation. 

Normalizing to the same /!# value (here choosen to be zero for simplicity) and 
comparing at the same initial energy density, the above mentioned 'strangeness 
enhancement' dissapears. 

The A s factor for particle collisions is within the error compatible with nucleus 
nucleus collisions as seen in figures |8| and [|. 

3 The increase is of course large when the threshold for strangeness production in p+p 
reactions is approached, while subthreshold strangeness production can still occur in nuclear 
collisions, but this limiting case is not the main point of the present discussion. 
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6) Enfin: what can we learn from strangeness production about the QCD phase 
transition ? 

First, figure |8| exhibits a very distinct feature of A s as a function of ef X s is in- 
creasing with initial energy density increasing from 0.14 to 1 GeV/fm 3 . After 
this €i value, X s saturates at a limiting value of A u m = 0.365 ± 0.033 ± 0.07 
I- 

Second, this behaviour is followed universally by all collisions studied. 
Third, A s is following closely the dependence of the temperature at the 
chemical freeze out of hadrons (as shown in figure |T0|). 

Assuming that local equilibrium conditions are reached in all collisions studied, 
and therefore that temperature is defined, we interpret the dramatic change 
of the temperature and of X s at chemical freeze-out near e« ~ 1 GeV/fm 3 , 
as critical behaviour: the onset of saturation occurs when the critical energy 
density for the QCD phase transition is reached. 

The T as well as X s would continue to rise with ej if no phase transition occurs. 

A bias in this interpretation is the implicit assumption that T does not cool 
down without noticing the transition due to some peculiar expansion dynamics 
and its y/s dependence. However in this case the cooling mechanism should 
generate the universal behaviour seen in figures ||] and |l^ by coincidence. 
Furthermore, since pressure and expansion characteristics have been measured 
to be different in different A+A collisions (e.g. there are strong flow phe- 
nomena pointing to high initial pressure in A+A collisions which are more 
pronounced at RHIC than in SPS, and not seen in p+p collisions) it seems 
difficult to explain a universally flat temperature curve exhibited by all p+p 
and A+A data, by a non universal expansion dynamics. 

A phenomenon analogous to the one seen in figure [K] is the following: 
We fill a box with water and look for the water-vapour phase transition with- 
out tools to detect vapour. Each time the transition to vapour (=QGP) occurs 
we thus wait until the vapour condensates back to water (=hadron gas), in 
order to measure its temperature. 

We make a plot of the water temperature as a function of the applied heat, 
and it looks as figure |l0j, namely it rises and saturates at the value of < 100° 
Celsius. 

Adding salt to water and repeating the experiment would result in different 
critical values rising with salinity. 

The baryochemical potential is like salt for hadronic systems. To achieve mea- 
suring one single curve one has to use the same salinity, as we do in figures 
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10] and H As a result, the border of the QCD phase transition can be drawn 
and the critical energy density can be extracted selfconsistently from the data, 
independent of any model predicting where the boundary must be. 



The dependence of the temperature on the baryon and strangeness chemical 
potentials is shown in figures [1] and |. 

If all these points would lie exactly on the critical border, their extrapolation 
to equivalent states at zero potentials would lead to one single critical temper- 
ature. 

If all points would not lie on the critical border but well below, their extrapo- 
lation to the equivalent states with \ib zero would give scattered temperatures 
below the critical one, without showing signs of saturation. 



One could counterargue that the saturation is caused by the freeze out con- 
ditions and the critical temperature is never reached in the initial state, (for 
example if there is no T cr n nearby). 

In this case however the freeze out conditions should change with the collision 
energy, because the cross sections of particle interactions and the mean free 
path changes too with \fs. Again saturation is then coincidental. 
A second argument against such coincidence is the following: If the freeze out 
conditions do imply a universal temperature at chemical freeze out then there 
should be no variation in the chemical freeze out temperature at zero chem- 
ical potentials. That is the plot of figure [TO should be completely flat for all 
thermal systems. 

A third argument is the similarity of the critical parameters that we extract 
in the above described picture and in [Q with the expected critical parameters 
from QCD (T mt ~ 200MeV, e crit ~ lGeV/fm 3 |§ and from other calcula- 
tions too e.g. from lattice simulations |TJ]). Even if the estimation of the initial 
energy density [^] has a large inherent error, e.g. in the assumption of a forma- 
tion time of 1 fm/c, still the error hardly amounts to a factor 7 (the difference 
between of Tevatron and the €i=l GeV//m 3 ). Also our limiting tempera- 
ture determination is more reliable than the initial energy density estimate and 
yet Tu m and e cr u agree reasonably with each other and with the QCD predic- 
tions. Finally we conclude that Bjorken's estimate |38| turns out to be correct. 



7) It is important to note, that while studing the \ s or the temperature de- 
pendence on the initial energy density using their values at non zero chemical 
potential (figure one can not extract the critical energy density neither the 
limiting temperature from this study. 

In particular the values one can extract in this case, simply reflect the fact 
that we deal with a critical border which is a 2-dimensional surface. However 
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more data are needed in order to observe this shift. 



8) All temperatures achieved with hadrons, are below the critical temperature 
for sure, however they can approach the border. We can measure only the 
approach to the critical surface from the hadronic side, study the way the 
critical values are approached and extract the critical parameters from critical 
behaviour. 



We study the way the systems approach the transition in figures [11], [12], [13 



In figures [14] and [Lg it is seen that the energy density over the particle density 
and the entropy density over the energy density show an onset of a critical 
behaviour at T ~ 150 MeV, while the entropy density over the particle density 
shown in figure |I5| saturates above T ~ 150 MeV. 

In figure [19] we try to fit the function / = (5 [ln{ 1 ./(l - T/T crit ))f (logarith- 
mic critical behaviour) to the data which go through the QGP phase transition 
according to figures |ll] and p (that is the points with T > 145), and extract 
the critical temperature T^ it and the critical exponent a^ lt from the fit. We 
find a value for the critical temperature of T^- t = 218 ± 70 MeV, and an expo- 
nent a fit = 0.54 ± 0.47 with \ 2 /DOF = 0.059/3. This behaviour can not be 
studied with the precision of a study of the neighbourhood of the Curie point 
in ferromagnets |3jJ . 

This value for T^ it is in agreement with the QCD expectations for the critical 
temperature of = 194 ± 18 MeV §. 



6 Conclusions 

The starting point of this paper is the extraction of thermodynamic parameters 
describing the final state of Au+Au collisions at 2 and 4 GeV per nucleon and 
of Pb+Pb collisions at 40 GeV per nucleon. We extrapolate these parameters 
to zero chemical potentials along an isentropic path and study the strangeness 
suppression factor A s (A s = as a function of the energy density reached 
early in each collision (initial energy density e^). 

We then arrive at the following conclusions: 

1) The so called 'strangeness suppression' puzzle, namely the decrease of the 
K/tt ratio (or equivalently of X s ) with a/s increasing from its value in Pb+Pb 
collisions at 40 A GeV, is explained as reflecting the varying chemical poten- 
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tials of the heavy ion systems. 

2) Several other experimental observations have the same origin: 

• the increase of the double ratio K/tt (A+A/p+p) with decreasing ^/s. 

• the flatter behaviour of the K/tt ratio as a function of ^/s when extracted 
at midrapidity. 

• the difference in the y/s dependence of K + /tt + and K~ /tt~ ratios. 

• the enhancement seen in strange particles in central p+A collisions as com- 
pared to p+p collisions at the same a/s. 

3) The X s value of systems with nonzero baryochemical potential is found to 
approach its limiting value at zero as defined by the pp and e + e~ colliding 
systems. We estimate that the limiting X s value (as well as an approximately 
net baryon free midrapidity region) will be reached by nuclear collisions at the 
initial energy density of ~ 8-9 GeV/fm 3 , (corresponding approximately to ^/s 
~ 3-8 TeV per nucleon+nucleon pair) probably at the LHC. 

4) Strangeness is not significantly increased in nucleus nucleus collisions as 
compared to elementary particle collisions, if they are compared 1) at the 
same (zero) chemical potential and 2) at the same initial energy density. 

5) However, X s is found to significantly increase in all systems which reach 
higher than ~ 1 GeV/fm 3 , as compared to all systems below. Strangeness 
is found to follow closely the temperature, rising until ~ 1 GeV/fm 3 and 
saturating along the border of the QCD phase transition, namely above 1 
GeV/fm 3 . This allows us to extract in a model independent way the critical 
parameters of the QCD phase transition from the data in particular e cr iu C ai = 
1 ± 0.3 GeV/fm 3 as well as the limiting T and X s values H. 

6) Having established in figures |^, and [TU] the critical initial energy density 
for the QCD phase transition of e crit ~ 1 ± 0.3GeV/ fm 3 , and determined the 
systems which go through the QGP phase, we study the way thermodynamic 
parameters approach the transition point. 

We find that the systems with > lGeV/ fm 3 , approach TA* = 218 ± 70 
MeV with a (logarithmic) critical exponent cr lt = 0.54 ± 0.47. 
More data from SPS, RHIC and LHC at several y^i will serve to narrow down 
the approach to the transition point. 
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Figure 1: The temperature as a function of the baryochemical potential for 
several nucleus+nucleus, hadron+hadron and lepton+lepton collisions. We 
demand for the fits confidence level > 10%. 
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Figure 2: The X s factor as a function of the baryochemical potential for several 
nucleus+nucleus, hadron+hadron and lepton+lepton collisions. We demand 
for the fits confidence level > 10%. 
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Figure 3: The baryochemical potential as a function of the strangeness chem- 
ical potential for several nucleus+nucleus collisions. We demand for the fits 
confidence level > 10%. 
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Figure 4: The temperature as a function of the strangeness chemical potential 
for several nucleus+nucleus hadron+hadron and lepton+lepton collisions. We 
demand for the fits confidence level > 10%. 
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Figure 5: The A s factor as a function of the strangeness chemical potential 
for several nucleus+nucleus hadron+hadron and lepton+lepton collisions. We 
demand for the fits confidence level > 10%. 
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Figure 6: (< A > + < K + K >)/n ratio as a function of F=f(A/i) in A+A 
and p+p collisions f28j, compared with a model from reference fl3l| . 
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Figure 7: The A s factor as a function of the initial energy density for several 
nucleus+nucleus, hadron+hadron and lepton+lepton collisions. We demand 
for the fits confidence level > 10%. The lines a and j3 show X s at nonzero 
while the line 7 show A s at zero 
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Figure 8: The \ s factor (in logarithmic scale) extrapolated to zero fugacities 
along an isentropic path, as a function of the initial energy density for several 
nucleus+nucleus, hadron+hadron and lepton+lepton collisions. We demand 
for the fits confidence level > 10%. 
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Figure 9: The X s factor (in linear scale) extrapolated to zero fugacities along 
an isentropic path, as a function of the initial energy density for several nu- 
cleus+nucleus, hadron+hadron and lepton+lepton collisions. We demand for 
the fits confidence level > 10%. 
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Figure 10: The temperature extrapolated to zero fugacities along an isentropic 
path, as a function of the initial energy density for several nucleus+nucleus, 
hadron+hadron and lepton+lepton collisions. We demand for the fits confi- 
dence level > 10%. 



27 




200 

T (MeV) 



Figure 11: The energy density p e as a function of the temperature extrap- 
olated to zero fugacities along an isentropic path for many nucleus+nucleus, 
hadron+hadron and lepton+lepton collisions. We demand for the thermal 
model fits confidence level > 10%. 
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Figure 12: The entropy density p s as a function of the temperature extrap- 
olated to zero fugacities along an isentropic path for many nucleus+nucleus, 
hadron+hadron and lepton+lepton collisions. We demand for the thermal 
model fits confidence level > 10%. 
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Figure 13: The density p n as a function of the temperature extrapo- 
lated to zero fugacities along an isentropic path for many nucleus+nucleus, 
hadron+hadron and lepton+lepton collisions. We demand for the thermal 
model fits confidence level > 10%. 
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Figure 14: The ratio of the energy density to the density as a function of the 
temperature extrapolated to zero fugacities along an isentropic path for many 
nucleus+nucleus, hadron+hadron and lepton+lepton collisions. We demand 
for the thermal model fits confidence level > 10%. 
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Figure 15: The ratio of the entropy density to the density as a function of the 
temperature extrapolated to zero fugacities along an isentropic path for many 
nucleus+nucleus, hadron+hadron and lepton+lepton collisions. We demand 
for the thermal model fits confidence level > 10%. 
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Figure 16: The ratio of the entropy density to the energy density as a function 
of the temperature extrapolated to zero fugacities along an isentropic path 
for many nucleus+nucleus, hadron+hadron and lepton+lepton collisions. We 
demand for the thermal model fits confidence level > 10%. 
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Figure 17: The ratio of the energy density to the temperature to the 4th power 
as a function of the temperature. The temperature is extrapolated to zero 
fugacities along an isentropic path for many nucleus+nucleus, hadron+hadron 
and lepton+lepton collisions. We demand for the thermal model fits confidence 
level > 10%. 
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Figure 18: The ratio of the entropy density to the temperature to the 
3th power as a function of the temperature. The temperature is extrapo- 
lated to zero fugacities along an isentropic path for many nucleus+nucleus, 
hadron+hadron and lepton+lepton collisions. We demand for the thermal 
model fits confidence level > 10%. 
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Figure 19: The X s factor as a function of the temperature both extrapo- 
lated to zero fugacities along an isentropic path for many nucleus+nucleus, 
hadron+hadron and lepton+lepton collisions. We demand for the thermal 
model fits confidence level > 10%. The lines shown are fits of the function 
/ = pin(l — T/T crit ) a below and above T=140 MeV. The fit of this func- 
tion to the data above T=140 MeV gives T crit = 218 ± 70 MeV, exponent 
a = 0.54 ± 0.47 and x 2 /DOF = 0.059/3. 
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